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Abstract 

This paper assesses the hedge effectiveness of an index-based longevity swap and a 
longevity cap. Although swaps are a natural instrument for hedging longevity risk, 
derivatives with non-linear pay-offs, such as longevity caps, also provide downside 
protection. A tractable stochastic mortality model with age dependent drift and 
volatility is developed and analytical formulae for prices of these longevity deriva¬ 
tives are derived. Hedge effectiveness is considered for a hypothetical life annuity 
portfolio. The hedging of the life annuity portfolio is comprehensively assessed for 
a range of assumptions for the market price of longevity risk, the term to maturity 
of the hedging instruments, as well as the size of the underlying annuity portfolio. 
The model is calibrated using Australian mortality data. The results provide a com¬ 
prehensive analysis of longevity hedging, highlighting the risk management benehts 
and costs of linear and nonlinear payoff structures. 


Key words: longevity risk management; longevity swaps; longevity options; hedge 
effectiveness 
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1 Introduction 


Securing a comfortable living after retirement is fundamental to the majority of the working pop¬ 
ulation around the world. A major risk in retirement, however, is the possibility that retirement 
savings will be outlived. Products that provide guaranteed lifetime income, such as life annuities, 
need to be offered in a cost effective way while maintaining the long run solvency of the provider. 
Annuity providers and pension funds need to manage the systematic mortality rislPI. associated 
with random changes in the underlying mortality intensity, in a life annuity or pension portfolio. 
Systematic mortality risk cannot be diversihed away with increasing portfolio size, while idiosyn¬ 
cratic mortality risk, representing the randomness of deaths in a portfolio with hxed mortality 
intensity, is diversihable. 


Reinsurance has been important in managing longevity risk for annuity and pension providers. 
However, there ar e concerns t h at rein surers have a limited risk appetite and are reluctant to take 
this “toxic” risk ( Blake et ahl (200^)). In fact, even if they were willing to accept the risk, the 
reinsurance sector is not deep enough to absorb the vast scale of longevity risk currently under¬ 
taken by annuity providers and pension fundsO] The sheer size of capital markets and an almost 
zero correlation between hnancial and demographic risks, suggests that they will increasingly take 
a role in the risk management of longevity risk. 


The hrst generatio n of capital market s olutio n s for longevi t y risk, in the form of mor t ality and 


longevity bonds ( iBlake and BurrowsI (1200 ll) . iBlake et ahl 
gained limited success. 


f 2006al) and iBauer et ahl f 2010lj lP1. 


The se cond generation involving forwards and swaps have attracted increasing interest (IBlake et al 
( 2ni3lH . Index-based instruments aim to mitigate systematic mortality risk, and hav e the poten 


tial to be less costly and are designed to allow trading as standardised contracts (IBlake et al. 
( 2ni3l H. Unlike the bespoke or customized hedging instrum ents such as re i nsura nce, they do not 


cover idiosyncratic mortality risk and give rise to basis risk (iLi and Hardvl (2011)). Since idiosyn¬ 
cratic mortality risk is reduced for larger portfolios, portfolio size is an important factor that 
determines the hedge effectiveness of index-based instruments. 


Longevity derivatives with a linear payoff, including q- forward s and S- forwards, have as an under¬ 
lying the mortality and the survival rate, respect ively (jLLMAl (j2010a|H . Their hedge effectiveness 
has been considered in iNgai and SherrisI (1201 ih who study the effectiveness of static hedging 
of longevity risk in different annuity portfolios. They consider a range of longevity-linked in¬ 
struments including q-forwards, longevity bonds and longevity swaps as hedging i nstruments to 
mitiga te longevity risk and demonstrate their benehts in reducing longevity risk. iLi and Hardv 
(I 2 OIII) also consider hedging longevity risk with a portfolio of q-forwards. They highlight basis 
risk as one of the obstacles in the development of an index-based longevity market. 


Longevi ty derivatives with a nonlin ear payoff structure have not received a great deal of attention 
to date. iBover and Stento'ft ( 20131) eval uate European and American type survivor options using 
simulations and IWang and Yang) ( 20131) propose and price survivor floors under an extension of 
the Lee-Carter model. These authors do not consider the hedge effectiveness of longevity options 
and longevity swaps as hedging instruments. 


^ From an annuity provider’s perspective, longevity risk modelling can lead to a (stochastically) over- or 
underestimation of survival probabilities for all annuitants. For this reason longevity risk is also referred 
to as the systematic mortality risk. 

^ It is estimated that pension assets for the 13 largest major pension markets have reached nearly 30 
trillions in 2012 (Global Pension Assets Study 2013, Towers Watson). 

^ Of particular interest is an attempt to issue the EIB longevity bond by the European Investment Bank 
(EIB) in 2004, which was underwritten by BNP Paribas. This bond was not well received by investors 


and cou ld not generate enough demand to be launched due to its deficiencies, as outlined in Blake et al. 
(l2(f06al L 
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Although dynamic hedging has been considered, because of the lack of liq uid markets in longevity 
risk, static hedging remains the only realistic option for annuity providers. ICairnsI (120111) considers 
q-forwards and a discrete-time delta hedging strategy, and compares it with static hedging. The 
lack of analytical formulas for pricing q-forwards and its derivatives, known as “Greeks”, can 
be a significant problem in assessing hedge effectiveness since simulations within simulations 
are required f or dynamic hedging strategies. The importance of tractable models has also been 
emphasis ed inlLuciano et al.l ( 2012 1 who also consider dynamic hedging for longevity and interest 
rate risk. Hari et al.l (1200811 apply a generalised two-factor Lee-Carter model to investigate the 
impact of longevity risk on the solvency of pension annuities. 


This paper provides pricing analysis of longevity derivatives, as well as their hedge effectiveness. 
We consider static hedging. A longevity swap and a cap are chosen as linear and nonlinear 
products to compare and assess index-based capital market products management of longevity 
risk management. The model used for this analysis is a continuous time model for mortality with 
age based drift and volatility, allowing tractable analytical formulae for pricing and hedging. The 
analysis is based on a hypothetical life annuity portfolio subject to longevity risk. The paper 
considers the hedging of longevity risk using a longevity swap and a longevity cap, a portfolio of 
S-forwards and longevity caplets respectively, based on a range of different underlying assumptions 
for the market price of longevity risk, the term to maturity of hedging instruments, as well as the 
size of the underlying annuity portfolio. 


The paper is organised as follows. Section 2 specifies the two-factor Gaussian mortality model, and 
its parameters are estimated using Australian males mortality data. Section 3 analyses longevity 
derivatives, in particular, a longevity swap and a cap, from a pricing perspective. Explicit pricing 
formulas are derived under the proposed two-factor Gaussian mortality model. Section 4 examines 
various hedging features and hedge effectiveness of a longevity swap and a cap on a hypothetical 
life annuity portfolio exposed to longevity risk. Section 5 summarises the results and provides 
concluding remarks. 


2 Mortality Model 


Let V 'Ht,P) be a filtered probability space where P is the real world probability 

measure. The subfiltration Qt contains information about the dynamics of the mortality intensity 
while death times of individuals are captured by l-if If is assumed that the interest rate r is 
constant where 5(0, t) = denotes the price of a f-year zero coupon bond, and our focus is on 
the modelling of stochastic mortality. 


2.1 Model Specification 


For the purpose of financial risk management applications one requires stochastic mortality model 
that is tractable, and is able to capture well the mortality dynamics for different ages. We work 
under the affine mortality intensity framework and assume the mortality intensity to be Gaus¬ 
sian such that analytical prices can be derived fo r longevity options , as described in Section 3. 
Gauss ian mortality models have been considered in Bauer et al.l(l2010h andiBlackburn and Sherris 
( 20131) within the forward mortality framework. Luciano and Vigna (2008) su ggest Gaussian mor - 


tality where the intensity follows the Ornstein-Uhlenbeck process. In addition, IJevtic et al.l (l2013l) 


consider a continuous time cohort model where the underlying mortality dynamics is Gaussian. 


We consider a two-factor Gaussian mortality model for the mortality intensity process pb^+tif) of 
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a cohort aged x at time t = ( 0 : 


dpx{t) = dYfit) + dY2{t), 

( 2 . 1 ) 

where 


dYfit) = aiYfit) dt + ai dWi{t) 

( 2 . 2 ) 

dY2{t) = {ax -\- fi) Y2{t) dt + ae^^ dW2{t) 

(2,3) 


and dWidW 2 = pdt. The first factor Ti(^) is a general trend for the intensity process that is 
common to all ages. The second factor Y 2 {t) depends on the initial age through the drift and 
the volatility termsThe initial values Ti(0) and 12 ( 0 ) of the factors are denoted by yi and y^, 
respectively. The model is tractable and for a specific choice of the par ameters (when a = 7 = 0) 
has been applied to short rate modelling in iBrigo and Mercnrio (2007). 


Proposition 2.1 Under the two-factor Gaussian mortality model (Eq. (2.1) - (2.3)j, the {T — t)- 
year expected survival probability of a person aged x + t at time t, conditional on filtration J^t, is 
given by 


s,+,((,r) =' e ! (e-Xt'-W*') = 


(2.4) 


where, using a 2 = ax and a 2 = ere'’'*, 


t ai{T-t) _ IN / 02(T-i) _ IN 

0(t,T) = - ’-Yfif) + - ’-Y 2 {t) and 


a\ 


02 


^ a? 


k=l “fe 




,2ofc(7’—i) 
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2pcri(j2 (^ ^ 0 _ X g“2(r 0 _ x ^ g(“i+«2)(T 0 _ x 

I 1 t T 


(y.\(y.2 


(y.\ 


0(2 


Oi + O 2 


(2.5) 


( 2 . 6 ) 


are the mean and the variance of the integral fj pbxiy)dv, which is Gaussian distributed, respec¬ 
tively. 

We will use the fact that the integral fj dv is Gaussian with known mean and variance to 
derive analytical pricing formulas for longevity options in Section 3. 

Proof. Solving Eq. (2.2) to obtain an integral form of Yi(t), we have 


Yfiu) du = dlTi(n)d 

The first term in Eq. (2.7) can be simplihed to 


(2.7) 




( g " i (^-0 - 1 ) 


CKi 


Yfit). 


^ For simplicity of notation we replace yx+t{t) by yx{t). 

^ We can in fact replace x by a; + t in Eq. (2.3). Using x -\- t will take into account the empirical 
observation that the volatility of mortality tends to increase along with age x -\-t (Figures 1 and 2). 
However, for a Gaussian process the intensity will have a non-negligible probability of reaching negative 
value when the volatility from the second factor becomes very high, which occurs for example 

when X -\-t > 100 (given 7 > 0). Using x instead of x +1 will also make the result in Section 3 easy to 
interpret. For these reasons we assume that the second factor Y 2 {t) depends on the initial age x only. 
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For the second term, we have 


cTi e-^^^dWi{v)du = dWi{v)du (— 

— r du fe"i“ r e-^^^dWi{v)) - — ( T e-"i’'dhFiO 

0^1 Jt \ Jt J (y.\ Jt \Jt 

r e-'^^^dWi{u) - — -IdWiiu), 

Jt ai Jt c^^ Jt 


O-l 

_ *^1 ^aiT 
0^1 J t 


(J^ P'^ 0*1 p^ 

— / e“i“e-“i“dhFi(M) = — / 

Jt 0C\ Jt 

where stochastic integration by parts is applied in the second equality. 


To obtain an integral representation for Y 2 {t), we follow the same steps as above, replacing Yi{t) 
by Y 2 {t) in Eq. (2.7). It is then straightforward to notice that 


.T 

fix{u) du = J Yi{u) + Y 2 {u) du 


( 2 . 8 ) 


is a Gaussian random variable with mean 0(t,T) (Eq. (2.5)) and variance r(t, T) (Eq. (2.6)). 
Equation (2.4) is obtained by applying the moment generating function of a Gaussian random 
variable. □ 


2.2 Parameter Estimation 


The discretised process, where the intensity is assumed to be const ant over each int e ger a ge 
and calendar year, is approximated by the central death rates m{x,t) ( Wills and SherrisI f 201lh ). 
Figure 1 displays Australian male central death rates m{x, t) for years t = 1970,1971,..., 2008 
and ages x = 60, 61,..., 95. Figure 2 shows the difference of the central death rates Am{x, t) = 
m{x + l,f + 1) — m{x,t). The variability of Am{x,t) is evidently increasing with increasing age 
X, which leads to the anticipation that 7 > 0. Furthermore, for a hxed age x, there is a slight 
improvement in central death rates for more recent years, compared to the past. 



Fig. 1. Australian male central death rates m{x, t) where t = 1970,1971,... , 2008 and x = 60,61,... , 95. 

The parameters {ai, a, 7 , pj, which det e rmine the volatility of the intensity process, are estimated 
as described below. As in iJevtic et al.l (1201311 . we aim to estimate parameters using the method 
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Fig. 2. Difference of the central death rates Am{x,t) = m{x + l,t + 1) — m{x,t) where 
t = 1970,1971, • • • , 2007 and x = 60,61, • • • , 94. 

of least squares, thus, calibrating the model to the mortality surface. However, we take advantage 
of the fact that a Gaussian model is employed where the variance of the model can be calculated 
explicitly and thus, we capture the diffusion part of the process by matching the variance of the 
model to mortality data. Specifically, the implemented procedure is as specified below: 

(1) Using empirical data for ages x = 60, 65,..., 90 we evaluate the sample variance of Am{x, t) 
across time, denoted by Var(Am3;). 

(2) The model variance Va.T{Afi^) for age x is given by 

Var(A/rJ = Var(criAlUi + ae^^AlUs) 

= (ai + 2aiape^^ + At. (2.9) 

Since the difference between the death rates is computed in yearly terms, we set At = 1. 

(3) The parameters {(Ti,cr, 7,p} are then estimated by fitting the model variance Vax^Ap^) to 
the sample variance Var(Ama;) for ages x = 60, 65,..., 90 using least squares estimation, 
that is, by minimising 

90 

(Var(A/i3.|(Ji,cr,7,p) - Var(Am3.))^ (2.10) 

x=60,65... 

with respect to the parameters {cxi, a, 7, p}. 

The remaining parameters {ai, a, P,yi,y 2 ^,yl^} are then estimated as described belo'vCll: 

(1) From the central death rates, we obtain empirical survival curves for cohorts aged 65 and 75 
in 2008. The survival curve is obtained by setting 

T 

Fa,(0,T) = ]^(1 -m(x + n - 1 , 0 )) ( 2 . 11 ) 

V=1 


® We calibrate the model for ages 65 and 75 simultaneously to obtain reasonable values for a and /3 
since the drift of the second factor ¥ 2 ( 1 ) is age-dependent. 
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where m{x, t) is the central death rate of an x years old at time tPI 
(2) The parameters {ai, a, /9, |/i, yl^} are then estimated by htting the survival curves (S'a;(0, T)) 
of the model to the empirical survival curves using least squares estimation, that is, by min¬ 
imising 

T 

E (2.12) 

x=65,75 j=l 

where Tes = 31 and T75 = 21, with respect to the parameters {ai, «,/9, |/i, 

The estimated parameters are reported in Table 1. Since 7 > 0 we observe that the volatility of 
the process is higher for older (initial) age x. 

Table 1 

Estimated model parameters. 

CJl (T 7 p Ol 

0.0022465 0.0000002 0.129832 -0.795875 0.0017508 

a 13 yi yf 

0.0000615 0.120931 0.0021277 0.0084923 0.0294695 


Age X = 65 Age x = 75 




Fig. 3. Percentiles of the simulated intensity processes and y 75 {t) for Australian males aged 65 

(upper left panel) and 75 (upper right panel) in 2008, with their corresponding survival probabilities 
(the mean and the 99% confidence bands) for a 65 years old (lower left panel) and 75 years old (lower 
right panel). 

The upper panel of Figure 3 shows the percentiles of the simulated mortality intensity for ages 
65 and 75 in the left and the right panel, respectively. One observes that the volatility of the 
mortality intensity is higher for a 75 year old compared to a 65 year old. Corresponding survival 
probabilities are displayed in the lower panel of Figure 3, together with the 99% conhdence bands 

^ Here t = 0 represents calendar year 2008 and we approximate the 1-year survival probability 

Q-m{x+v-l,0) Py y _ 
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computed pointwise. As it is pronounced from the figures, the two-factor Gaussian model specihed 
above, despite its simplicity, produces reasonable mortality dynamics for ages 65 and 75. 


3 Analytical Pricing of Longevity Derivatives 


We consider longevity derivatives with different payoff structures including longevity swaps, 
longevity caps and longevity floors. Closed form expressions for prices of these longevity deriva¬ 
tives are derived under the assumption of the two-factor Gaussian mortality model introduced 
in Section 2. These instruments are written on survival probabilities and their properties are 
analysed from a pricing perspective. 


3.1 Risk-Adjusted Measure 


For the purpose of no-arbitrage valuation, we require the dynamics of the factors Yi{t) and Y 2 {t) to 
be written under a risk-adjusted measure[E| To preserve the tractability of the model, we assume 


that the processes Wi{t) and W 2 {t) with dynamics 

dWiit) = dWiit) (3.1) 

dW2{t) = Xae^^ Y2{t) dt + dW2{t) (3.2) 

are standard Brownian motions under a risk-adjusted measure Q. In Eq. (3.2) A represents the 
market price of longevity riskj^] Under Q we can write the factor dynamics as follows: 

dYiit) = aiYi{t) dt -|- ui dWi{t) (3.3) 

dY 2 {t) = (ax + /3 — Xae'^^) 12(f) dt (J 2 dW 2 (t). (3.4) 

The corresponding risk-adjusted survival probability is given by 


where a 2 = ax(3 replaced by (ax + /S — Xae'^^) in the expressions for 0(f, T) and T(t, T), see 
Eq. (2.5) and Eq. (2.6), respectively. 


Since a liquid longevity market is yet to be developed, we aim to determine a reasonable value 
for A based on the longevity b ond announc e d by BNP Paribas an d European Investmen t Ban k 
(EIB) i n 2004 as propose d in ICairns et ^ (l2006l) and applied in iMevricke and SherrisI f 2014il . 
see also IWills and SherrisI (2011). The BNP/EIB longevity bond is a 25-year bond with coupon 
payments linked to a survivor index based on the realised mortality ratesjZ^ The price of the 
longevity bond is given by 


25 


r(0) = B(0, T) 

T=1 


(e-Io 


dv 


(3.6) 


° Since the longevity market is still in its development stage and hence, incomplete, we assume a risk- 
adjusted measure exists but is not unique. 

® For simplicity, we assume that there is no risk adjustment for the first factor Yi and A is age- 
independent. 

The issue price was determined by BNP Paribas using anticipated cash flows based on the 2002-based 
mortality projections provided by the UK Government Actuary’s Department. 















where <5 is a spread, or an average risk premium per annumF^. and the T-year projected survival 
rate is assumed to be the T-year survival probability for the Australian males cohort aged 65 as 
modelled in Section 2, see Eq. (2.4). S ince the BNP/EIB bond is priced based on a yield of 20 
basis points below standard EIB rates f Cairns et al.l (1200611 1. we have the spread of <5 = 0.00203 


Under a risk-adjusted measure 


the price of the longevity bond corresponds to 


25 


J2B{0,T)g, 

T=1 


,Q(A) 

0 


Jo t^x{v)dv 


(3.7) 


Fixing the interest rate to r = 4%, we hnd a model-dependent A, such that the risk-adjusted bond 
price U‘®^^^(0) matches the market bond price U(0) as close as possible. For example, for A = 8.5 
we have U(0) = 11.9045 and U '^l'^HO) = 11.9068. For more details on the above procedure refer 
to Mevricke and SherrisI ( 20141) . In the following we assume that the risk-adjusted measure 
determined by a unique value of A. 


IS 



Fig. 4. Risk-adjusted survival probability with respect to different market price of longevity risk A. 

Figure 4 shows the risk-adjusted survival probabilities for Australian males aged 65 with respect 
to different values of the market price of longevity risk A. As one observes from the hgure, a larger 
(positive) value of A leads to an improvement in survival probability, while a smaller values of A 
indicate a decline in survival probability under the risk-adjusted measure Q. 


3.2 Longevity Swaps 


A longevity swap involves counterparties swapping hxed payments for payments linked to the 
number of survivors in a re ference popu lation in a given time period, and can be thought of as a 
po rtfolio o f S-forw ards, see iDow^ ( 2003f) . An S-forward, or ‘survivor’ forward has been developed 
by ILLMAI ( 2010b|) . Longevity swaps can be regarded as a stream of S-forwards with different 


The spread 6 de pends on the term of the bond and the initial age of the cohort being tracked 
( Cairns et al.l ( 20061 11. and 6 is related to but distinct from A, the market price of longevity risk. 

The reference cohort for the BNP/EIB longevity bond is the England and Wales males aged 65 in 
2003. Since the longevity derivatives market is under-developed in Australia, we assume that the same 
spread of (i = 0.002 (as in the UK) is applicable to the Australian males cohort aged 65 in 2008. Note 
however that sensitivity analyses will be performed in Section 4. 
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maturity dates. One of the advantages of using S-forwards is that there is no initial capital 
requirement at the inception of the contract and cash flows occur only at maturity. 

Consider an annuity provider who has an obligation to pay an amount dependent on the number 
of survivors, and hence, survival probability of a cohort at time T. If longevity risk is present, the 
survival probability is stochastic. In order to protect himself from a larger-than-expected survival 
probability, the provider can enter into an S-forward contract paying a hxed amount K G (0,1) 
and receiving an amount equal to the realised survival probability exp { — /q dv} at time T. 

In doing so, the survival probability that the provider is exposed to is certain, and corresponds 
to some hxed value K. If the contract is priced in such a way that there is no upfront cost at the 
inception, it must hold that 


B{0,T)E^ 



r 

Jo 


fix{v) dv 



(3.8) 


under the risk-adjusted measure Q. Thus, the hxed amount can be identihed to be the risk- 
adjusted survival probability, that is. 


K{T) = B? 



(3.9) 


Assuming that there is a positive market price of longevity risk, the longevity risk hedger who pays 
the hxed leg a nd receives the bo ating leg bears the cost for entering an S-forwardFollowing 
terminology in Bifhs et 1 20141 ). the amount K{T) = ^3,(0, T) can be referred to as the swap 


rate of an S-forward with maturity T. In general, the mark-to-market price process F(t) of an 
S-forward with hxed leg K (not necessarily K{T) as in Eq. (3.9)) is given by 


F{t) = B{t,T)Ff 
= B(t,T)FP 

= B{t, T) (^3,(0, t) Sx+t{t, T) - k) (3.10) 

for f G [0,T]. The quantity 

is the realised survival probability, or the survivor index for the 

The term Sx+t{t,T) that appears in Eq. (3.10) has a natural interpretation. Given in¬ 

formation Fq at time f = 0, this term becomes S'a;(0,T), which is the risk-adjusted survival 
probability. As time moves on and more information Ft-, with t G (0,T), is revealed, the term 
Sx{0,t) Sx+tit,T) is a product of the realised survival probability of the hrst t years, and the 
risk-adjusted survival probability in the next (T — t) years. At maturity T, this product becomes 
the realised survival probability up to time T. In order words, one can think of ^3,(0, t) Sx+t(t, T) 
as the T-year risk-adjusted survival probability with information known up to time t. 


(3.11) 

cohort, which is observable given 


The price process F{t) in Eq. (3.10) depends on the swap rate Sx+tit,T) of an S-forward written 
on the same cohort that is now aged {x -|- t) at time t, with time to maturity {T — t). If a 
liquid longevity market was developed, the swap rate Sx+tit,T) could be obtained from market 
data. As S'a;(0,t) is observable at time t, the mark-to-market price process of an S-forward could 
be considered model-independent. However, since a longevity market is still in its development 
stage, market swap rates are not available and a model-based risk-adjusted survival probability 

The risk-adjusted survival probability will be larger than the “best estimate” P-survival probability if 
a positive market price of longevity risk is demanded, see Figure 4. 
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Sx+t{t,T) has to be used instead. An analytical formula for the mark-to-market price of an S- 
forward can be obtained if the risk-adjusted survival probability is expressed in a closed-form, 
which can be performed, for example, under the two-factor Gaussian mortality model. 

Since a longevity swap is constructed as a portfolio of S-forwards, the price of a longevity swap 
is simply the sum of the individual S-forward prices. 


3.3 Longevity Caps 


A longevity cap, which is a portfolio of longevity caplets, provides a similar hedge to a longevity 
swap but is an option-type instrument. Consider again a scenario described in Section 3.2 where 
an annuity provider aims to hedge against larger-than-expected T-year survival probability of a 
particular cohort. Alternatively to hedging with an S-forward, the provider can enter into a long 
position of a longevity caplet with payoff at time T corresponding to 

max I (^e-ro| (3.12) 

where K G (0,1) is the strike priceOE If the realised survival probability is larger than K, the 
hedger receives an amount (^exp {—/q^/ i3;(n) dn} — from the longevity caplet. This payment 
can be regarded as a compensation for the increased payments that the provider has to make 
in the annuity portfolio, due to the larger-than-expected survival probability. There is no cash 
outflow if the realised survival probability is smaller than or equal to K. In other words, the 
longevity caplet allows the provider to “cap” its longevity exposure at K with no downside risk. 
Since a longevity caplet has a non-negative payoff, it comes at a cost. The price of a longevity 
caplet 

Ci{t-,T,K) = B{t,T)E^ ('313) 

under the two-factor Gaussian mortality model is obtained in the following Proposition. 

Proposition 3.1 Under the two-factor Gaussian mortality model (Eg. {2.1)-Eg. (2.3)j the price 
at time t of a longevity caplet Ci{t;T, K), with maturity T and strike K, is given by 

Clit- T, K) = St St Bit, T) ^ [\]Tit,T) - d) - KBit, T)$ (-d) (3.14) 


where St = S^i^S) is the realised survival probability observable at time t, St = Sx+tit,T) is the 
risk-adjusted survival probability in the next (T —f) years, d = ,} fin {iP/(^f^t)} -|- |f(t, T)^ 

and $(•) denotes the cumulative distribution function of a standard Gaussian random variable. 


Proof. Under the risk-adjusted measure Q, we have, from Proposition (2.1), that 

L j%xiv)dv ~ iV(-0(t,T),f(f,T)). 

Using the simplified notation 0 = 0(t, T), P = P(f, T) we can write 
C^it■, T, K) = Bit, T)E^ [iSt - K)+) 

roo 1 _l( Gh® V / _ 

= Bit,T) I ^=e HvAJ 


(3.15) 


\/27rP 


Ste^-K] de 


The payoff of a longevity caplet i s simil ar to the payoff of the option embedded in the principal-at-risk 
bond described in Biffis and Blakd (2014). 
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Vr 


1 


di-K r ^ . 


Equation (3.14) follows using properties of <|)(-) and noticing that St = ”, that is, 0 = 

if-ln^*. □ 


Similar to an S-forward, the price of a longevity caplet depends on the product term ^a;(0, t) Sx+tit, T). 
In particular, a longevity caplet is said to be out-of-the-money \i K > Sx{0,t) Sx+tit, T); at-the- 
money a K = ^3,(0, t) Sx+tit, T)] and in-the-money if < ^3,(0, t) ^3;+t(t, T). Eq. (3.14), is verihed 
using Monte Carlo simulation summarised in Table 2, where we set r = 4%, A = 8.5 and t = 0. 
Other parameters are as specihed in Table 1. 

Table 2 

Pricing longevity caplet Ci{Q',T,K) by the formula (Eq. (3.14)) and by Monte Carlo simulation of 
Eq. (3.13); [ , ] denotes the 95% confidence interval. 


(T, K) 

Exact 

M.C. Simulation 

(10, 0.6) 

0.15632 

0.15644 [0.15631, 0.15656] 

(10, 0.7) 

0.08929 

0.08941 [0.08928, 0.08954] 

(10, 0.8) 

0.02261 

0.02262 [0.02250, 0.02275] 

(20, 0.3) 

0.08373 

0.08388 [0.08371, 0.08406] 

(20, 0.4) 

0.03890 

0.03897 [0.03879, 0.03914] 

(20, 0.5) 

0.00525 

0.00530 [0.00522, 0.00539] 


Following the result of Proposition 3.1, the two-factor Gaussian mortality model leads to the price 
of a longevity caplet that is a function of the following variables: 

• realised survival probability ^3,(0, f) of the hrst t years; 

• risk-adjusted survival probability Sx+tit, T) in the next T — t years; 

• interest rate r; 

• strike price K; 

• time to maturity (T — f); and 

• standard deviation (f, T), which is a function of the time to maturity and the model param¬ 
eters. 

Since the quantity exp | — /(J' /ix('y) dnj is log-normally distributed under the two-factor Gaussian 
mortality model, Eq. (3.14) resembles the Black-Scholes formula for option pricing where the 
underlying stock price follows a geometric Brownian motion. In our setup, the stock price at time 
t is replaced by the T-year risk-adjusted survival probability ^3,(0, t) Sx+tit, T) with information 
available up to time t. While the stock is traded and can be modelled directly using market data, 
the underlying of a longevity caplet is the survival probability which is not tradable but can be 
determined as an output from the dynamics of mortality intensity. As a result, the role of the stock 
price volatility in the Black-Scholes formula is played by the standard deviation of the integral 
of the mortality intensity ^xiv)dv. Since the integral ^xiv) dv captures the whole history 
of the mortality intensity ^xit) from t to T under Q, one can interpret the standard deviation 

(t, T) as the volatility of the risk-adjusted aggregated longevity risk of a cohort aged x + t aX 
time t, for the period from t to T. 

The left panel of Figure 5 shows caplet prices for a cohort aged x = 65, using parameters as 
specihed in Table 1, as a function of time to maturity T and strike K. We set r = 0.04, A = 8.5 
and t = 0 such that S'3,(0, 0) = 1. A lower strike price indicates that the buyer of a caplet is willing 
to pay more to secure a better protection against a larger-than-expected survival probability. On 
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Fig. 5. Caplet price as a function of (left panel) T and K and (right panel) A where K = 0.4 and T = 20. 

the other hand, when the time to maturity T is increasing, the underlying survival probability is 
likely to take smaller values, which leads to a higher probability for the caplet to become out-of- 
the-money at maturity for a hxed iC, see Eq. (3.12). Consequently, for a hxed K the caplet price 
decreases with increasing T. 

The right panel of Figure 5 illustrates the effect of the market price of longevity risk A on the 
caplet price. The price of a caplet increases with increasing A. As shown in Figure 4, a larger 
value of A will lead to an improvement in survival probability under Q. Thus, a higher caplet 
price is observed since the underlying survival probability is larger (on average) under Q when A 
increases, see Eq. (3.13). 

Since longevity cap is constructed as a portfolio of longevity caplets, it can be priced as a sum of 
individual caplet prices, see also Section 4.1.2. 


4 Managing Longevity Risk in a Hypothetical Life Annuity Portfolio 


Hedging features of a longevity swap and cap are examined for a hypothetical life annuity portfolio 
subject to longevity risk. Factors considered include the market price of longevity risk, the term 
to maturity of hedging instruments and the size of the underlying annuity portfolio. 


4-1 Setup 


We consider a hypothetical life annuity portfolio that consists of a cohort aged a; = 65. The size 
of the portfolio that corresponds to the number of policyholders, is denoted by n. The underlying 
mortality intensity for the cohort follows the two-factor Gaussian mortality model described in 
Section 2, and the model parameters are specihed in Table 1. We assume that there is no loading 
for the annuity policy and expenses are not included. 

Further, we assume a single premium, whole life annuity of $1 per year payable in arrears condi¬ 
tional on the survival of the annuitant to the payment dates. The fair value, or the premium, of 
the annuity evaluated at t = 0 is given by 


a,= ^R(0,T)A,(0,T) (4.1) 

r=i 

where r = 4% and u = 110 is the maximum age allowed in the mortality model. The life annuity 
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provider, thus, receives a total premium, denoted by A, for the whole portfolio corresponding to 
the sum of individual premiums: 

A = nttx- (4-2) 

This is the present value of the asset held by the annuity provider at t = 0. Since the promised 
annuity cashflows depend on the death times of annuitants in the portfolio, the present value of 
the liability is subject to randomness caused by the stochastic dynamics of the mortality intensity. 
The present value of the liability for each policyholder, denoted by L^, is determined by the death 
time Tk of the policyholder, and is given by 


[rfc J 

Lk=Y.B{0,T) (4.3) 

T=1 

for a simulated r^, with |_gj denoting the next smaller integer of a real number q. The present 
value of the liability L for the whole portfolio is obtained as a sum of individual liabilities: 

n 

L = J2Lt. (4.4) 

k=l 

The algorithm for simulating death times of annuitants, which requires a single simulated path 
for the mortality intensity of the cohort, is summarised in Appendix A. The discounted surplus 
distribution (i^no) of an unhedged annuity portfolio is obtained by setting 

Dno = A - L. (4.5) 

The impact of longevity risk is captured by simulating the discounted surplus distribution where 
each sample is determined by the realised mortality intensity of a cohort. Since traditional pricing 
and risk management of life annuity relies on diversihcation effect, or the law of large numbers, 
we consider the discounted surplus distribution per policy 

Dno/n. (4.6) 

Figure 6 shows the discounted surplus distribution per policy without longevity risk (i.e. when 
setting (Ji = cr = 0) with different portfolio sizes, varying from n = 2000 to 8000. As expected, 
the mean of the distribution is centred around zero as there is no loading assumed in the pricing 
algorithm, while the standard deviation diminishes as the number of policies increases. 

In the following we consider a longevity swap and a cap as hedging instruments. These are 
index-based instruments where the payoffs depend on the survivor index, or the realised survival 
probability (Eq. (3.11)), which is in turn determined by the realised mortality intensity. We do 
not consider basis risM^ but due to a hnite portfolio size, the actual proportion of survivors, 
, where Nt denotes the number of deaths experienced by a cohort during the period [0 ,t], 
will be in general similar, but not identical, to the survivor index (Appendix A). As a result, the 
static hedge will be able to reduce systematic mortality risk, whereas the idiosyncratic mortality 
risk component will be retained by the annuity provider. 


4-1.1 A Swap-Hedged Annuity Portfolio 

For an annuity portfolio hedged by an index-based longevity swap, payments from the swap 

'e-Io >^Av)dv _ (4.7) 


If basis risk is present, we need to distinguish between the mortality intensity for the population (pj!) 
and mortality intensity for the cohort (p^) underlying the annuity portfolio, see Biffis et al. ( 20141 ). 
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Discounted Surplus Distribution Per Policy 



Fig. 6. Discounted surplus distribution per policy without longevity risk with different portfolio size (n). 

at time T G ,T} depend on the realised mortality intensity, where T denotes the term to 
maturity of the longevity swap. The number of policyholders n acts as the notional amount of 
the swap contract so that the quantity nexp{— dv} represents the number of survivors 

implied by the realised mortality intensity at time T. We hx the strike of a swap to the risk- 
adjusted survival probability, that is, 

K{T) = S^{0,T)= (4.8) 

such that the price of a swap is zero at f = 0, see Section 3.2. The discounted surplus distribution 
of a swap-hedged annuity portfolio can be expressed as 

D = A - L + F 14 9j 

■^swap ±j ^ ± swap 


where 

Fswap = nj2 B{0,T) - S,{0,T)) (4.10) 

T=1 ^ ^ 

is the (random) discounted cashflow coming from a long position in the longevity swap. The 
discounted surplus distribution per policy of a swap-hedged annuity portfolio is determined by 

hi swap/ 

4-1.2 A Cap-Hedged Annuity Portfolio 

For an annuity portfolio hedged by an index-based longevity cap, the cashflows 

n max | - K{T)^ , o| (4.11) 

at T G {1,...,T} are payments from a long position in the longevity cap. We set 

K{T) = S^{0,T) = (4.12) 
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such that the strike for a longevity caplet is the “best estimated” survival probability given 
The discounted surplus distribution of a cap-hedged annuity portfolio is given by 
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-Dcap — A — L + Fcap — Ccap (4-13) 

where 

Fcap = n ^ 5(0,T) - S',(0,T)) ,o| (4.14) 

is the (random) discounted cashflow from holding the longevity cap and 

f 

C,,p = nJ2Ci{0;T,S,{0,T)) (4.15) 

T=1 

is the price of the longevity cap. The discounted surplus distribution per policy of a cap-hedged 
annuity portfolio is given by Dca.p/n. 


4-2 Results 


Hedging results are summarised by means of summary statistics that include mean, standard 
deviation (std. dev.), skewness, as well as Value-at-Risk (VaR) and Expected Shortfall (ES) of 
the discounted surplus distribution per policy of an unhedged, a swap-hedged and a cap-hedged 
annuity portfolio. Skewness is included since the payoff of a longevity cap is nonlinear and the 
resulting distribution of a cap-hedged annuity portfolio is not symmetric. VaR is defined as the 
g-quantile of the discounted surplus distribution per policy. ES is dehned as the expected loss 
of the discounted surplus distribution per policy given the loss is at or below the g-quantile. We 
£x g = 0.01 so that the confidence interval for VaR and ES corresponds to 99%. We use 5,000 
simulations to obtain the distribution for the discounted surplus. Hedge effectiveness is examined 
with respect to (w.r.t.) different assumptions underlying the market price of longevity risk (A), 
the term to maturity of hedging instruments (T) and the portfolio size (n). Parameters for the 
base case are as specified in Table 3. 

Table 3 

Parameters for the base case. _ 


A 

T (years) 

n 

8.5 

30 

4000 


4-2.1 Hedging Features w.r.t. Market Price of Longevity Risk 

The market price of longevity risk A is one of the factors that determines prices of longevity 
derivatives and life annuity policies. Since payoffs of a longevity swap, a cap and a life annuity are 
contingent on the same underlying mortality intensity of a cohort, all these products are priced 
using the same A. Figure 7 and Table 4 illustrate the effect of changing A on the distributions 
of an unhedged, a swap-hedged and a cap-hedged annuity portfolio. The degree of longevity risk 
can be quantihed by the standard deviation, the VaR and the ES of the distributions. We observe 
that increasing A leads to the shift of the distribution to the right, resulting in a higher average 
surplus. On the other hand, changing A has no impact on the standard deviation and the skewness 
of the distribution. 


For a longevity swap, the risk-adjusted survival probability is used as a strike price so that the price 
of a longevity swap is zero at inception. In contrast, a longevity cap has non-zero price and Sx{0,T) is 
the most natural choice for a strike. 
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Fig. 7. Effect of the market price of longevity risk A on the discounted surplus distribution per policy. 


Table 4 


Hedging features of a longevity swap and cap w.r.t. market price of longevity risk A. 



JUean 

Std.dev. 

Skewness 

VaRo.99 

ESo.99 


A = U 

JNo hedge 

-0.0076 

0.3592 

-0.2804 

-0.9202 

-1.1027 

Swap-hedged 

-0.0089 

0.0718 

-0.1919 

-0.1840 

-0.2231 

Cap-hedged 

-0.0086 

0.2054 

1.0855 

-0.3193 

-0.3515 


A = 4.5 

JNo hedge 

0.1520 

0.3592 

-0.2804 

-0.7606 

-0.9431 

Swap-hedged 

0.0048 

0.0718 

-0.1919 

-0.1703 

-0.2094 

Cap-hedged 

0.0682 

0.2054 

1.0855 

-0.2425 

-0.2746 


A = 8.5 

JNo hedge 

[12978 

0:8592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

0.0204 

0.0718 

-0.1919 

-0.1547 

-0.1938 

Cap-hedged 

0.1205 

0.2054 

1.0855 

-0.1903 

-0.2224 


A = 12.5 

Mo hedge 

0.4475 

0.3592 

-0.2804 

-0.4650 

-0.6476 

Swap-hedged 

0.0398 

0.0718 

-0.1919 

-0.1354 

-0.1744 

Cap-hedged 

0.1619 

0.2054 

1.0855 

-0.1489 

-0.1810 


For an unhedged annuity portfolio, a higher A leads to higher premium for the life annuity 
policy since the annuity price is determined by the risk-adjusted survival probability Sx(0,T), 
see Eq. (4.1). In other words, an increase in the annuity price compensates the provider for the 
longevity risk undertaken when selling life annuity policies. There is also a trade-off between risk 
premium and affordability. Setting a higher premium will clearly improve the risk and return of an 
annuity business, it might, however, reduce the interest of potential po licyholders. An empiric al 


relationship between implied longevity and annuity prices is studied in IChigodaev et al.l (1201411 . 


When life annuity portfolio is hedged using a longevity swap, the standard deviation and the 
absolute values of the VaR and the ES reduce substantially. The higher return obtained by 
charging a larger market price of longevity risk in life annuity policies is offset by an increased 
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price paid implicitly in the swap contract (since Sx{0,T) > Sx{0,T) in Eq. (4.10)). It turns out 
that as A increases an extra return earned in the annuity portfolio and the higher implicit cost of 
the longevity swap nearly offset each other out on average. The net effect is that a swap-hedged 
annuity portfolio remains to a great extent unaffected by the assumption on A, leading only to a 
very minor increase in the mean of the distribution. 

For a cap-hedged annuity portfolio, the discounted surplus distribution is positively skewed since 
a longevity cap allows an annuity provider to get exposure to the upside potential when policy¬ 
holders live shorter than expected. Compared to an unhedged portfolio, the standard deviation 
and the absolute values of the VaR and the ES are also reduced but the reduction is smaller com¬ 
pared to a swap-hedged portfolio. When A increases, we observe that the mean of the distribution 
for a cap-hedged portfolio increases faster than for a swap-hedged portfolio but slower than for 
an unhedged portfolio. It can be explained by noticing that when the survival probability of a 
cohort is overestimated, that is, when annuitants turn out to live shorter than expected, hold¬ 
ing a longevity cap has no effect (besides paying the price of a cap for longevity protection at 
the inception of the contract) while there is a cash outflow when holding a longevity swap, see 
Eq. (4.10) and Eq. (4.14). 

In the longevity risk literature, the VaR and the ES are of a particular importance as they are 
the main factors determinin g the capital reserve when dealing with exposure to longevity risk 
( Mevricke and SherrisI fl2014i) b As shown in Table 4, the difference between a swap-hedged and 
a cap-hedged portfolio in terms of the VaR and the ES becomes smaller when A increases. In 
fact, for A > 17.5, a longevity cap becomes more effective in reducing the tail risk of an annuity 
portfolio compared to a longevity swapjZI] This result suggests that a longevity cap, besides being 
able to capture the upside potential, can be a more effective hedging instrument than a longevity 
swap in terms of reducing the VaR and the ES when the demanded market price of longevity risk 
A is large. 


4-2.2 Hedging Features w.r.t. Term to Maturity 
Table 5 


Hedging features of a longevity swap and cap w.r.t. term to maturity T. 



JVlean 

Std.dev. 

Skewness 

VaKo.99 

ESo,99 


r = 10 Years 

JNo hedge 

012978 

013592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

0.2820 

0.2911 

-0.3871 

-0.5707 

-0.7490 

Cap-hedged 

0.2893 

0.2989 

-0.2661 

-0.5801 

-0.7592 


T = 20 Years 

JNo hedge 

012978 

013592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

0.1740 

0.1794 

-0.7507 

-0.3656 

-0.5061 

Cap-hedged 

0.2234 

0.2310 

0.2006 

-0.3870 

-0.5259 


T = 30 Years 

JNo hedge 

0.2978 

0.3592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

0.0204 

0.0718 

-0.1919 

-0.1547 

-0.1938 

Cap-hedged 

0.1205 

0.2054 

1.0855 

-0.1903 

-0.2224 


T = 40 Years 

JNo hedge 

0.2978 

0.3592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

-0.0091 

0.0668 

0.0277 

-0.1616 

-0.1869 

Cap-hedged 

0.0984 

0.1999 

1.1527 

-0.1909 

-0.2131 


Table 5 and Figure 8 summarize hedging results with respect to the term to maturity of hedging 
instruments. Due to the long-term nature of the contracts, the hedges are ineffective for T < 10 
years and the standard deviations are reduced only by around 17 — 19% for both instruments. 
The lower left panel of Figure 3 shows that there is little randomness around the realised sur¬ 
vival probability for the first few years for a cohort aged 65, and consequently the hedges are 


Given A = 17.5, the VaR and the ES for a swap-hedged portfolio are —0.1051 and —0.1441 respectively. 
For a cap-hedged portfolio they become —0.1038 and —0.1360, respectively. 
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Fig. 8. Effect of the term to maturity T of the hedging instruments on the discounted surplus distribution 
per policy. 

insignificant when T is short. 

The difference in hedge effectiveness between T = 30 and T = 40 for both instrnments is also 
insignihcant. In fact, the longevity risk nnderlying the annnity portfolio becomes small after 30 
years since the majority of annnitants has already deceased before reaching the age of 95. In onr 
model setnp the chance for a 65 years old to live np to 95 is aronnd 6% (Fignre 4 with A = 0) 
and, hence, only aronnd 4000 x 6% = 240 policies will still be in-force after 30 years. Much of 
the risk left is attributed to idiosyncratic mortality risk, and hedging longevity risk for a small 
portfolio using index-based instruments is of limited use. 

For a swap-hedged portfolio, the standard deviation is reduced signihcantly when T > 20 years. 
The mean surplus, on the other hand, drops to nearly zero since there is a higher cost implied for 
the hedge with increasing number of S-forwards involved to form the swap as T increases. 

Similar hedging features with respect to T are observed for a longevity cap. However, the skewness 
of the distribution of a cap-hedged portfolio increases with increasing T. It can be explained by 
noticing that while a longevity cap is able to capture the upside potential regardless of T, it 
provides a better longevity risk protection when T is larger. As a result, the distribution of a 
cap-hedged portfolio becomes more asymmetric when T increases. 


4-2.3 Hedging Features w.r.t. Portfolio Size 

Table 6 and Figure 9 demonstrate hedging features of a longevity swap and a cap with changing 
portfolio size n. We observe a decrease in standard deviation, as well as the VaR and the ES (in 
absolute terms) when portfolio size increases. Compared to an unhedged portfolio, the reduction 
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Fig. 9. Effect of the portfolio size n on the discounted surplus distribution per policy. 


Table 6 


Hedging features of a longevity swap and cap w.r.t. different portfolio size (n). 



iVlean 

Std.dev. 

Skewness 

VaKo.99 

RSo.99 


n = 2UUU 

JNo hedge 

[12973 

013046 

-0.2662 

-0.6360 

-0.8107 

Swap-hedged 

0.0200 

0.0990 

-0.1615 

-0.2120 

-0.2653 

Cap-hedged 

0.1200 

0.2160 

0.9220 

-0.2432 

-0.2944 


n = 4000 

JNo hedge 

012978 

013592 

-0.2804 

-0.6148 

-0.7973 

Swap-hedged 

0.0204 

0.0718 

-0.1919 

-0.1547 

-0.1938 

Cap-hedged 

0.1205 

0.2054 

1.0855 

-0.1903 

-0.2224 


n = 0000 

Mo hedge 

012977 

013566 

-0.2786 

-0.6363 

-0.8001 

Swap-hedged 

0.0204 

0.0594 

-0.3346 

-0.1259 

-0.1660 

Cap-hedged 

0.1204 

0.2016 

1.1519 

-0.1639 

-0.2051 


n = 8000 

Mo hedge 

012982 

013554 

-0.2920 

-0.6060 

-0.7876 

Swap-hedged 

0.0209 

0.0536 

-0.5056 

-0.1190 

-0.1595 

Cap-hedged 

0.1209 

0.1992 

1.1616 

-0.1598 

-0.1991 


in the standard deviation and the risk measures is larger for a swap-hedged portfolio, compared 
to a cap-hedged portfolio. Recall that idiosyncratic mortality risk becomes signihcant when n 
is small. We quantify the effect of the portfolio size on hedge effectiveness by introducing the 
measure of longevity risk reduction R, dehned in terms of the variance of the discounted surplus 
per policy, that is. 


Var(R)*) 
Var(R)) ’ 


(4.16) 


where Var(Z)*) and Var(Z)) represent the variances of the discounted surplus distribution per 
policy for a hedged and an unhedged annuity portfolio, respectively. The results are reported in 
Table 7. 
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Table 7 

Longevity risk reduction R of a longevity swap and cap w.r.t. different portfolio size (n). 


n 

2000 

4000 

6000 

8000 

-^swap 

92.6% 

96.0% 

97.2% 

97.7% 

-Reap 

64.9% 

67.3% 

68 .0% 

68 .6% 


iLi and Hardvl f 201ll) consider hedging longevity risk using a portfolio of q-forwards and find the 
longevity ri sk reduction of 77.6% and 69.6% for portfolio size of 10,000 and 3,000, respectively. In 


contrast to iLi and Hardvl (120 llh . we do not consider basis risk and the result of using longevity 
swap as a hedging instrument leads to a greater risk reduction. Overall, our results indicate that 
hedge effectiveness for an index-based longevity swap and a cap diminishes with decreasing n 
since idiosyncratic mortality risk cannot be effectively diversihed away for a small portfolio size. 
Even though a longevity cap is less effective in reducing the variance, part of the dispersion is 
attributed to its ability of capturing the upside of the distribution when survival probability of 
a cohort is overestimated. From Table 6 we also observe that the distribution becomes more 
positively skewed for a cap-hedged portfolio when n increases, which is a consequence of having 
a larger exposure to longevity risk with increasing number of policyholders in the portfolio. 


5 Conclusion 


Life and pension annuities are the most important types of post-retirement products offered by 
annuity providers to help securing lifelong incomes for the rising number of retirees. While interest 
rate risk can be managed effectively in the hnancial markets, longevity risk is a major concern 
for annuity providers as there are only limited choices available to mitigate the long-term risk. 
Development of effective hnancial instruments for longevity risk in capital markets is arguably 
the best solution available. 


Two types of longevity derivatives, a longevity swap and a cap, are analysed in this paper from a 
pricing and hedging perspective. We apply a tractable Gaussian mortality model to capture the 
longevity risk, and derive explicit formulas for important quantities such as survival probabilities 
and prices of longevity derivatives. Hedge effectiveness and features of an index-based longevity 
swap and a cap used as hedging instruments are examined using a hypothetical life annuity 
portfolio exposed to longevity risk. 


Our results suggest that the market price of longevity risk A is a small contributor to hedge 
effectiveness of a longevity swap since a higher annuity price is partially offset by an increased 
cost of hedging when A is taken into account. It is shown that a longevity cap, while being able to 
capture the upside potential when survival probabilities are overestimated, can be more effective 
in reducing longevity tail risk compared to a longevity swap, provided that A is large enough. 
The term to maturity T is an important factor in determining hedge effectiveness. However, the 
difference in hedge effectiveness is only marginal when T increases from 30 to 40 years for an 
annuity portfolio consisting of a single cohort aged 65 initially. This is due to the fact that only 
a small number of policies will still be in-force after a long period of time (30 to 40 years), and 
index-based instruments turn out to be ineffective when idiosyncratic mortality risk becomes a 
larger contributor to the overall risk, compared to systematic mortality risk. The effect of the 
portfolio size n on hedge effectiveness is quantihed and compared with the result obtained in 
iLi and Hardvl ( 20111) where population basis risk is taken into account. In addition, we hnd that 
the skewness of the surplus distribution of a cap-hedged portfolio is sensitive to the term to 
maturity and the portfolio size, and, as a result, the difference between a longevity swap and a 
cap when used as hedging instruments becomes more pronounced for larger T and n. 


As discussed in iBifSs and Blakel ( 2014i) . developing a liquid longevity market requires reliable 
and well-designed hnancial instruments that can attract sufficient amount of interests from both 
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buyers and sellers. Besides of a longevity swap, which is so far a common longevity hedging 
choice for annuity providers, option-type instruments such as longevity caps can provide hedging 
features that linear products cannot offer. A longevity cap is shown to have alternative hedging 
properties compared to a swap, and this option-type instrument would also appeal to certain 
classes of investors interested in receiving premiums by selling a longevity insurance. Further 
research on the design of longevity-linked instruments from the perspectives of buyers and sellers 
would provide a further step towards the development of an active longevity market. 


A Appendix 


To simulate death times of annuitants, we notice that once a sample of the mortality intensity is 
obtained, the Cox process becomes an inhomogeneous Poisson process and t he hrst jump times, 
which are interpreted as death times, can be simulated as follows (see e.g. iBrigo and Mercuriol 

(j2003)): 

(1) Simulate the mortality intensity /ix(t) from t = 0 to t = u) — x. 

(2) Generate a standard exponential random variable For example, using an inverse transform 
method, we have — In (1 — u) where u ~ Uniform(0,1). 

(3) Set the death time r to be the smallest T such that ^ fix{s) ds. If ^ fix{s) ds 

then set T = u) — X. 

(4) Repeat step (2) and (3) to obtain another death time. 


The payoff of an index-based hedging instrument depends on the realised survival probability 
exp {—/q Hx{v)dv}. The payoff of a customised instrument, on the other hand, depends on the 
proportion of survivors, , underlying an annuity portfolio where the number of deaths, Nf, 
is obtained by counting the number of simulated death times that are smaller than t. Note that 


3“/o 


n — N+ 


n 


(A.l) 


and the accuracy of the approximation improves when n increases. 
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